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(SJ ■ ABSTRACT 

We study the effects of primordial non-Gaussianity on the large scale structure in 
the excursion set approach, accounting for correlations between steps of the random 
• walks in the smoothed initial density field. These correlations are induced by realistic 

\ smoothing filters (as opposed to a filter that is sharp in fc-space) , but have been ignored 

by many analyses to date. We present analytical arguments - building on existing 
results for Gaussian initial conditions - which suggest that the effect of the filter at 
large smoothing scales is remarkably simple, and is in fact identical to what happens 
in the Gaussian case: the non-Gaussian walks behave as if they were smooth and 
deterministic, or "completely correlated" . As a result, the first crossing distribution 
(which determines, e.g., halo abundances) follows from the single-scale statistics of 
■ the non-Gaussian density field - the so-called "cloud- in-cloud" problem does not exist 

for completely correlated walks. Also, the answer from single-scale statistics is simply 
£f) . one half that for sharp- k walks. 

' We explicitly test these arguments using Monte Carlo simulations of non-Gaussian 

, walks, showing that the resulting first crossing distributions, and in particular the 

^sO ■ factor 1/2 argument, are remarkably insensitive to variations in the power spectrum 

and the defining non-Gaussian process. We also use our Monte Carlo walks to test 
some of the existing prescriptions for the non-Gaussian first crossing distribution. 
OO ' Since the factor 1/2 holds for both Gaussian and non-Gaussian initial conditions, it 

' provides a theoretical motivation (the first, to our knowledge) for the common practice 

of analytically prescribing a ratio of non-Gaussian to Gaussian halo abundances. 
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1 INTRODUCTION tial scales, and the study of the former therefore promises 

_ . . _ . ... _ . . ,,,„,. , , , to complement the information gained on primordial NG 

the detection of primordial non-Gaussiamty NG) (Malda- „ , , n . , , ., , , ,. „ nl , , 

, . from the latter, bmcc detailed studies ot collapse and struc- 
cena 2003; Acquaviva et al. 2003) would serve as a power- „ , . . , . , .. 

' ture formation involve expensive numerical simulations, it 

ml discriminator between models for seeding curvature in- . . . n . n n . n 

, . , it 1S interesting to explore the extent to which one can gain 



analytical insight into the problem. 



homogeneities in the early universe (Babich, Creminelli & 
Zaldarriaga 2004; LoVerde et al. 2008; Sartoris et al. 2010; 

Sefusatti 2009; see Desjacques & Seljak 2010 for a recent re- The excursion set approach (Epstein 1983; Bond et al. 

view). Traces of primordial NG can potentially be found not 1991; Lacey & Cole 1993) has gone a long way in helping to 

only in the cosmic microwave backround (CMB) radiation build such insight. This approach simplifies the problem of 

(Creminelli et al. 2006; Komatsu et al. 2010), but also in the structure formation by separating its highly non-linear, dy- 

late time large scale structure of the universe (Slosar et al. namical aspects from its statistical ones. The key idea which 

2008). Both the abundance of collapsed objects (halos) and the approach hinges on is the realisation that matter does 

correlations in their spatial distribution are sensitive to the not move very far while halos form. This allows the problem 

presence of NG in the initial conditions (see, e.g., Dalai et of the late time distribution of matter, to be mapped onto 

al. 2008; Matarrese & Verde 2008; Afshordi & Tolley 2008). its initial conditions. Whether or not a structure forms on 

Collapsed objects and the CMB probe very different spa- a given length scale is decided by asking whether or not the 

initial density field, smoothed on that scale, lies above a cer- 
tain threshold (Press & Schechter 1974). The density field 

* E-mail: musso@ictp.it performs a random walk as the smoothing scale is changed, 

f E-mail: aparanja@ictp.it The excursion set ansatz relates the statistics of these ran- 
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dom walks to those of large scale structure (e.g., halo abun- 
dances, bias, merger rates, formation times, etc.), while the 
dynamics enters through the choice of the threshold, which 
acts as an absorbing barrier for the walks. Most importantly 
for the present discussion, the mapping of the problem onto 
the initial conditions brings primordial NG into the fold of 
the excursion set approach. 

Halo abundances - the focus of this paper - in the ex- 
cursion set approach are prescribed in terms of the first 
crossing distribution (absorption rate) of a barrier, by the 
random walks performed by the smoothed density field. If 
dn/dm is the differential number density of halos with mass 
in (m, m + dm), then 

™^) dm=/(s)ds . (i) 
p dm 

where p is the mean matter density of the universe, s is the 
variance of the density field on scale R with m = AivpR 3 /3, 
and /(s) is the first crossing distribution. In the context 
of Gaussian initial conditions, it was recognised early on 
(Bond et al. 1991) that the first crossing distribution f(s) 
admits an analytical solution when the field is smoothed 
with a filter that is sharp in fc-space, while the problem is 
much more involved with a more realistic choice of filter, 
such as, e.g., a TopHat in real space. While Bond et al. 
presented numerical results (which we discuss below), it is 
worth noting that Peacock & Heavens (1990) had already 
presented an analytical approximation to the first crossing 
distribution of filtered walks, which described the numerical 
results of Bond et al. remarkably accurately. 

Based on Eq. {TJ, it certainly appears desirable to calcu- 
late a first cross distribution, rather than perform an analy- 
sis analogous to Press & Schechter (1974) by only calculating 
unconditional probabilities for the smoothed density con- 
trast to be above threshold (although see below). In the con- 
text of non-Gaussian initial conditions, however, the techni- 
cal hurdles involved have meant that several analyses have 
relied on extending the original Press-Schechter argument to 
the non-Gaussian case (Matarrese, Verde & Jimenez 2000; 
LoVerde et al. 2008; LoVerde & Smith 2011). Only recently 
have there been systematic efforts to tackle the problem 
using the full excursion set approach of Bond et al. (e.g., 
Lam & Sheth 2009; Maggiore & Riotto 2010; D'Amico et al. 
2011a), very few of which have explicitly addressed the issue 
of non-sharp-fc filtering (e.g., Maggiore & Riotto 2010, who 
work with the TopHat filter). 

One issue which has hampered the understanding of 
non- Gaussian halo abundances, is the fact that even for 
Gaussian initial conditions, both types of analyses (Press- 
Schechter and Bond et al.) are known not to describe well the 
results of TV-body simulations. In part this is because the ex- 
ponential tail of the halo mass function is not well described 
by the canonical threshold value S c — 1.686 from spherical 
collapse (Gunn & Gott 1972), a point we will return to later. 
The standard practice in the literature has therefore been 
to use Press-Schechter/sharp-fc analyses to prescribe a ratio 
Rng = /NG//gauss (e.g. LoVerde et al. 2008). As yet, there 
has been no convincing explanation of why such a prescrip- 
tion can be expected to work, although JV-body simulations 
with non-Gaussian initial conditions indicate that the pre- 
scription does work (Grossi et al. 2007; Pillepich, Porciani & 
Hahn 2008; Desjacques & Seljak 2010). And finally, it is also 



worth keeping in mind that equating the halo abundances 
on the left hand side of Eq. Jl} with the first crossing distri- 
bution on the right, is only an ansatz (see Paranjape, Lam 
& Sheth 2011a, henceforth PLS11, for a discussion). 

Overall then, the present situation as regards the non- 
Gaussian halo mass function is, in our opinion, rather 
clouded. There are several calculations of the function 
/ng(s) which appears on the right hand side of Eq. fT}. Ac- 
cepting that /ng (s) must be a first crossing distribution, we 
see that none of the calculations of /ng(s) have been tested 
using numerical first crossing distributions. Consequently, 
one cannot claim to understand which (if any) prescription 
for /ng(s) works well. The ratio prescription has been tested 
in TV-body simulations and seems to do well; however, as yet 
there is no theoretical understanding of why this must be 
so. There have been very few attempts to assess the impact 
of nontrivial filtering on the calculation of /ng(s), mainly 
based on a linearization of the problem (Maggiore & Riotto 
2010). However, since the expansion parameter involved is of 
order ~ 0.4, we would argue that the actual effect of filtering 
on non- Gaussian random walks is not yet clear. 

Our aim in this paper is to gain some insight into these 
issues. We will therefore focus on testing analytical pre- 
scriptions for non-Gaussian first crossing distributions using 
Monte Carlo simulations of random walks (as opposed to 
mass functions from iV-body simulations), along the lines 
initiated by PLS11 for the Gaussian case. We will demon- 
strate that the first crossing distribution of sharp-fc non- 
Gaussian walks on the relevant scales is rather insensitive to 
the exact nature of the defining non-Gaussian process, and 
depends, instead, mainly on the values of the leading order 
parameters (such as the skewness) which define the process. 
Our main result, however, concerns the effects of a nontriv- 
ial filter. We will present analytical arguments, supported by 
our numerics, to show that for small values of variance s the 
first crossing distribution of filtered non-Gaussian walks is 
simply a factor 1/2 times the corresponding distribution for 
sharp-fc walks. This effect is also remarkably robust against 
dramatic changes in, e.g., the power spectrum and type of 
non-Gaussianity. Since the same effect is known to be true 
for the case of Gaussian walks, as an immediate consequence 
we see that the ratio /ng / /gauss of first crossing distributions 
at small s is independent of the choice of filter. While this 
does not on its own explain the effect being seen in sim- 
ulations, since there is a big leap involved in going from 
first crossing distributions to mass functions, we believe it 
at least makes the ratio prescription theoretically plausible. 

The paper is organized as follows. In section 2 we 
present our analytical arguments, explaining the role of 
the filter. Section 3 deals with the numerical solution; we 
show results for different choices of power spectra and non- 
Gaussian processes. A final section summarizes and dis- 
cusses some avenues for further work. The Appendix con- 
tains some technical details regarding the structure of non- 
Gaussianity. 

2 ANALYTICAL ARGUMENTS: A FACTOR 

OF 1/2 

The problem of evaluating the first crossing distribution 
of random walks with correlated steps for Gaussian initial 
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conditions, was solved numerically by Bond et al. (1991). 
In the cosmological context, it is the linearly extrapolated, 
smoothed density field 8 which performs a random walk in 
steps of its variance 

s(R) = {5 2 n} = J dkk 2 ^W 2 (kR) , (2) 

where Ps(k) is the power spectrum of matter fluctuations, 
linearly evolved to present day, and W(kR) is a smooth- 
ing filter. While a filter that is sharp in fc-space leads to 
walks with uncorrelated steps, which can be treated ana- 
lytically, more realistic filters such as the TopHat (in real 
space) or Gaussian lead to nontrivial correlations between 
steps. Among other things, Bond et al. demonstrated that 
the high-mass or small s tail of the first crossing distribution 
/(s) is very well described, independently of the choice of fil- 
ter and power spectrum, by the Press-Schechter expression 

r 1 \ 1 <5c -&l/2s fn\ 

s/psOO = o^= e ' ( 3 ) 

^ V 27TS 

which is a factor 1/2 different from the analytical answer for 
walks with uncorrelated steps. The analytical approximation 
presented by Peacock & Heavens (1990) also reduces to this 
expression for small s, independently of the filter and power 
spectrurrQ. In both cases, the actual details of the filter and 
power spectrum only come into play when the relation be- 
tween s and the smoothing scale R is needed. 

This behaviour can be understood in terms of what 
PLS11 recently discussed as the "completely correlated" 
limit of the random walks. In this limit, the walks are de- 
terministic in the sense that each walk has a constant value 
of v = 5/y/s. In the plane of 6/y/s vs. s, each walk therefore 
crosses the barrier 5 = 5 C at most once. The "survival proba- 
bility", i.e. the fraction of walks that has not crossed at s, is 
therefore simply the fraction of walks for which v < 8 c /yfs, 
which is (1 + erf (5 c /v / 2s))/2 for Gaussian initial conditions. 
The derivative of this survival probability leads to the Press- 
Schechter expression ([3]) , which is therefore the first crossing 
distribution for walks with completely correlated steps and 
Gaussian initial conditions. 

In other words, the so-called cloud-in-cloud problem 
which Bond et al. solved for the sharp-fc case, does not ex- 
ist for completely correlated walks; such walks either cross 
the constant barrier exactly once, or not at all. For corre- 
lations induced by realistic filters (such as the TopHat or 
Gaussian), Eq. Q remains a good description at small s, 
since the walks have not had enough "time" to start feeling 
the effects of stochasticity; in the small s limit, the walks 
are still completely correlated. 

Now consider the case of non- Gaussian initial condi- 
tions. In what follows it will be useful to define the nor- 
malised equal-scale connected moments of the non-Gaussian 
field 8, €j = ( 5 J+2 } c /s" , the most important for us be- 
ing the three-point function given by ei, 

f ~ {53{S)) (4) 

61 = s 3/2 ■ W 

1 The choice of filter does affect how long the small- s regime lasts. 
As one can see from Fig. 9 of Bond et al. (see also Paranjapc, 
Lam & Sheth 2011a), the completely correlated answer is accurate 
at least up to ^/s ~ 0.5 <5 C for both TopHat and Gaussian filters. 



A recent analysis by Maggiore & Riotto (2010) has explored 
TopHat filter effects in the calculation of the first crossing 
distribution, using a path-integrals analysis. In contrast, we 
will be guided by the discussion above of the case of Gaus- 
sian initial conditions. In particular, for small s that entire 
argument goes through as is for the non-Gaussian case as 
well, apart from the fact that, at any given s, the fraction of 
walks Png(< <5c/%/*) with v < 8 c /yfs is no longer simply an 
error function, but receives non-Gaussian corrections. These 
corrections, however, can be calculated simply by knowing 
the single-pomt statistics of the non-Gaussian filtered field 8, 
even though this field in general has a complicated structure 
of unequal-scale correlations. This happens because, as far 
as the tail of the first crossing distribution is concerned, any 
memory induced by the non-Gaussianity is rendered super- 
fluous by these filter-induced strong correlations. We check 
this explicitly in the next section, by numerically calculating 
the first crossing distribution of non-Gaussian walks. 

To summarize, we have argued that in its small s tail 
the first crossing distribution of walks with non-Gaussian 
initial conditions and filter-induced correlated steps must 
be the derivative of the survival probability for walks with 
completely correlated steps: 

d d 
s/ng(s) = -s— P surv (s) -> -s— P NG (< 8 c /\fs) ■ (5) 

We emphasize that the final result is valid only for small s; 
the result for Gaussian initial conditions (Peacock & Heav- 
ens 1990; Bond et al. 1991) suggests that the shape of the dis- 
tribution at large s could be very different. This is not a ma- 
jor concern, since the effect of primordial non-Gaussianity 
is expected to be large only at small s, and consequently we 
need not worry about deviations from this simple result at 
large s. 

This final result is not new in itself; the single-point 
statistics of the non-Gaussian field have been the basis of 
the calculations of, e.g., Matarrese, Verde & Jimenez (2000) 
(henceforth MVJ00), LoVerde et al. (2008) and Lo Verde 
& Smith (2011). However, since those analyses ultimately 
aimed at obtaining a mass function rather than a first cross- 
ing distribution, the true statistical meaning of those calcu- 
lations has become clouded, in our opinion, and is clarified 
by our arguments above. 

It is worth comparing the structure of the sharp-fc non- 
Gaussian first crossing distribution with that of the uncon- 
ditional probability Png(< 8 c /i/s) in some detail, since this 
yields some insight into the ratio prescription mentioned 
in the Introduction. The calculation of Maggiore & Riotto 
(2010) shows that the effect of unequal scale correlations of 
the non-Gaussian field 8 in the first crossing distribution, 
are typically suppressed by powers of u^ 1 = \fs/8 c in the 
small s tail. As we explain below, if we therefore ignore these 
contributions altogether, then the derivative with respect to 
s of the probability Png(< S c /y/s), when expressed in terms 
of the equal-scale connected moments tj, will be exactly 1/2 
times the sharp-fc first crossing distribution, 

- S^Png(< 8 c /t/s) — i X s/ N G,sharp-fc(s) • (6) 

The reason is that in each case, sharp-fc as well as completely 
correlated, the non-Gaussianity appears as an exponentiated 
derivative operator ~ exp[— eic^/6 + • • •] acting on the cor- 
responding Gaussian first crossing distribution. This can be 
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seen, e.g., by comparing equation (A3) of Paranjape, Gor- 
don & Hotchkiss (2011) (henceforth PGH11) for the sharp-fc 
case, with equations (26) and (36) of MVJ00 for the com- 
pletely correlated case. 

The theoretical effort now lies in computing the ef- 
fect of this derivative operator, and different groups have 
approached the problem with different assumptions (see 
e.g. MVJ00, LoVerde et al. 2008, Maggiore & Riotto 2010, 
D'Amico et al. 2011a, PGH11). The bottom-line for us is the 
same in each case: since the Gaussian results for the sharp- 
fc and completely correlated cases differ by a factor of 1/2, 
consequently so do the non-Gaussian answers. An immediate 
consequence is that at small s, the ratio of non-Gaussian to 
Gaussian first crossing distributions /NG//gauas will be the 
same regardless of whether the walks were sharp-fc filtered 
or, more realistically, completely correlated. 

As with the case of Gaussian initial conditions, the rela- 
tion between s and smoothing scale R explicitly depends on 
the choice of filter and power spectrum. Additionally now, 
the moments tj will also in general be sensitive to these 
choices. What is universal though, is the functional form of 
the right hand side of Eq. © when expressed in terms of s 
and the ej. 

For later comparison, we note here that the completely 
correlated prediction derived from the PGH11 prescription^ 2 ] 
reads (with v = Jc/y^s) 

S/PGH,cc(s) 



2 ^27r(l + eii/) 

t\v — (1 + eiu) ln(l + t\u) 



x exp 



while the corresponding prediction for MVJ00 is 



(7) 



s/MVJ,cc(s) = - 



^(1 - ieii/) 



(8) 

in which we assumed ei to be constant, which is what we 
will use for our numerical results below. 



3 MONTE CARLO TESTS 

To verify that filter-induced correlations do in fact com- 
pletely correlate the walks at small s, in this subsection 
we numerically calculate the first crossing distributions for 
non-Gaussian random walks, comparing the effect of sharp-fc 
and Gaussian filters. The results are not expected to change 
when switching to the TopHat filter; we choose the Gaussian 
filter since it allows us to derive simple analytical expressions 
for some of the quantities we will need below. Along the way 
we will also demonstrate that, for the sharp-fc walks as well, 
the first crossing distribution is rather insensitive to the de- 
tails of the unequal scale correlations of the non-Gaussian 
field. 



3.1 Gaussian initial conditions 

To set the stage, let us begin by discussing the case of Gaus- 
sian initial conditions in some detail. The strategy for gen- 



2 I.e., 1/2 times their expression for the sharp-fc sf(s), the latter 
being their equation (A9) divided by 2. 
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Figure 1. Distribution of scale s at which a walk first crossed 
the constant barrier <5 = 8 C (histograms), for sharp-fc (solid) and 
Gaussian filtered (dashed) walks with Gaussian initial conditions. 
While the filtered walks assumed a power law power spectrum 
Pg(k) oc fc — 12 , the result is insensitive to the choice of P&(k). The 
histograms agree with the analytical predictions (smooth curves) 
of Bond et al. (1991) (solid) and Eq. (J3j (dashed), respectively. 



crating walks with Gaussian initial conditions is the same as 
that employed by Bond et al. (1991): the sharp-fc walks are 
constructed by accumulating independent Gaussian random 



numbers gi with a fixed variance As, 5 



(sharp — k) 



while the Gaussian filtered walks are obtained by applying 
the filter W(kR) = e~ < - kR) /2 to the same set of numbers to 
get dj Gau) = J2i gtWihRj). In the latter case, one also needs 
to know which values of kj and Rj to associate to the j-th 
step. Once a power spectrum is specified, this can be done by 
inverting the relations jAS = (2tt 2 )~ 1 J^ 3 dkk 2 P$(k) and 
jAS= (2-K 2 )- 1 J+°° dkk 2 P s (k)W 2 (kR j ). 

In this work, we will use a power law form P$(k) oc fc n 
for the power spectrum, rather than the ACDM form, in or- 
der to simplify the analysis. The results of Bond et al, and 
our previous analytical arguments, suggest that the final re- 
sult for the non-Gaussian first crossing distribution should 
be independent of the form of the power spectrum. We will 
also explicitly verify this expectation by displaying results 
for two very different values of n. As a sanity check, the his- 
tograms in Figure [T]show our numerical solution for the first 
crossing distribution of the constant barrier 8 = 5 C for walks 
with Gaussian initial conditions, for the sharp-fc and Gaus- 
sian filters respectively, with n = —1.2. The corresponding 
smooth curves show the analytical expressions of Bond et al. 
(1991) and Eq. ([3}, respectively, which differ by a factor 1/2. 
The agreement indicates that our code is working correctly. 
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3.2 Non-Gaussian initial conditions: model 
definitions 

To generate non-Gaussian walks, we must choose a type of 
non-Gaussianity. Ideally, we would pick as the initial density 
field a realization of a non-Gaussian process compatible with 
the bispectrum shape of some specific model of inflation, and 
smooth it with the filter of choice to obtain the walks. For 
example, in the so-called local model (Komatsu & Spergel 
2001; Lyth, Ungarelli & Wands 2003; Bartolo, Matarrese & 
Riotto 2004; Dvali, Gruzinov & Zaldarriaga 2004) the pri- 
mordial curvature perturbation £(x) is written in real space 
as £ ~ + /nlCgi where £g is a Gaussian random field. 
The density contrast S is then obtained from £ applying the 
Poisson equation and some matter transfer function. 

In practice, however, the nonlinear part of this 8 in- 
volves a convolution in Fourier space over the modes of a 
Gaussian field So, which is somewhat cumbersome to han- 
dle. Since our goal is mainly to demonstrate the effect of the 
filter (namely, the factor of 1/2 in the first crossing distri- 
bution), it is sufficient to adopt a simpler model. Our main 
requirement will be that our model must reproduce an equal- 
scale three point function e\ of our choice. In the local model, 
ei is weakly scale dependent; for our numerical walks, we will 
simply choose a constant value for t\. Additionally, we will 
also consider specific choices of unequal-scale behaviour for 
the three-point function. We discuss the unequal scale be- 
haviour of the three-point function of the local model, in 
Appendix Al. 

We will consider the class of integrated non- Gaussian 
processes defined by 



S(S) = 5 G {S) + 



ds 



8 G (S)S G (s) - (Sg(S)Sg(s)) 

tO) 



(9) 



where Sg is the Gaussian field smoothed with the filter of 
choice, the dot is a derivative with respect to the argument 
s, and P and 7 are functions that can be tuned to reproduce 
the desired behavior of the three-point function. As we show 
in Appendix A2, setting 



P(S) = eiVS 



s dc (6 G (S)S G (s)} 
< 7(«) 



(10) 



allows us to match exactly any value of ei, for any choice 
of filter. This expression for /3(S) ignores corrections from a 
term involving e 3 ; we have checked that, for the values of ei 
we consider in this work, including this correction does not 
affect any of our numerical results. For later use, note that 
for a sharp-fc filter, (Sg(S)Sg(s)) = 1, while for the Gaussian 
filter and a power law Ps(k) oc k", we have (Sg(S)Sg(s)) = 

2 (n+3)/2j 1 + / s /g\2/(n+3)]-(n+S)/2_ 

Additionally, setting 'y(s) — 4 v / s approximately repro- 
duces the unequal scale behavior of {8(Si)5(S2)8(Ss)) of 
the filtered local model, and makes the model above a good 
proxy to study non-Gaussianity of the local type (see Ap- 
pendix A2). We emphasize, however, that the details of the 
non-Gaussian process are actually expected to be irrelevant 
for the first crossing distribution. One could, e.g., pick the 
function f(s) so as to generate any other choice of unequal 
scale behaviour. Below, we will also display the results of 
simply setting 7(5) = 1. 

We will refer to the case 7(5) = 4t/s as model A. When 



generating the random walks for this model, we find it more 
convenient to analytically perform an integration by parts 
and bring the model definition to the forrrjj 



S(S) = S a (S) + 



Cl 



6(1 + a) 



Sg(S) 2 



+ 5 G (S) 



ds 

2s 3 / 2 



(11) 



where a = 1 for the sharp-fc filter, and for the Gaussian 
filter with Pj(fe) oc k n , a = (1.52, 1.31) for n = (-1.2, -2), 
which follows from using the expression for ( Sq(S)8g(s) ) 
mentioned earlier, and performing the resulting integrals. 

The second model we will consider (henceforth model 
B) is defined by setting 7(5) = 1 in Eq. ((9J and choosing j3 
appropriately, 



8(s) = S G (s) + 



Sg(s) 



(12) 



(see Appendix IA1 for details). This model significantly dif- 
fers from model A, and hence from the local model as well, 
in the structure of the unequal scale three point function 
( 5(si)8(s2)5{sz) ) . As we show in the Appendix, in the limit 
when si < S2 < S3, model B has ( 8(si)8(s2)8(sz) ) oc 
Si^/s2 whereas the local model has ( 8(sx)8{s%)8(s3) }i oc oc 

V / SlS2. 

Model B is also interesting because its single-point 
statistics can be calculated analytically (see, e.g., Eqs. 2- 
5 of MVJ00). We do this in Appendix B and show that the 
first crossing distribution for walks with completely corre- 
lated steps (the analog of Eqs. [3] [7]and[8]l is 



f f> 1 Sc 

SJNG, modB(S) = ~ — ?=PmodB 
2 <JS 



where, defining X(v, ei) = (ei/3)(i> + ei/6) 
1 



(13) 



PmodB^, Cl) = 



: exp 



-JLfl + A^er)) 



V27r(l + 2A(z/,ci)) 

x 2 cosh r^ v /l + 2A(i/,ei)^ . (14) 

We will use this below to test our claim that single-point 
statistics are sufficient to determine the first crossing distri- 
bution of filtered walks. 

3.3 Non-Gaussian walks: sharp-fc case 

Our emphasis here is on the factor 1/2 which arises due to 
the walks being completely correlated at small s. As such, it 
would be enough to simply display the ratio of the first cross- 
ing distributions for Gaussian filtered and sharp-fc walks, 

3 Comparing this form with the approximate description of the 
local model in Eq. (11) of Afshordi and Tolley (2008), we see 
that our model A is in fact quite close to the local model. 
Their term involving Q p G&mG can he compared with our term 
~ S(S) f dss~ 3 / 2 8(s): at least with Gaussian filtering, one has 
8(R) ~ Jd 3 kk 2 ^e- k2R2 ^ ~ (d/dR 2 )<S>(R), so that *(s) ~ 
J dsg(s)8(s) for some suitable choice of g(s). This is another way 
of understanding why, as discussed in the Appendix, the unequal 
scale behaviour of the three point function of model A is close to 
that of the local model. 
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Figure 2. Distribution of scale s at which a walk first crossed 
the constant barrier 8 = S c (histograms), for sharp-fc walks with 
non-Gaussian initial conditions defined by models A (solid) and 
B (dotted), for ei = 0.3 (red, upper) and ei = —0.3 (blue, lower). 
Note that the histograms for the different models are remark- 
ably alike. The smooth curves show the theoretical predictions of 
PGH11 (solid) and MVJ00 (dotted), with the red (upper) and 
blue (lower) curves corresponding to positive and negative ei, re- 
spectively, with | e 1 1 = 0.3 in each case. See text for a discussion. 



checking that it is consistent with being 1/2, and we do this 
below. We would like to do more, however, and also compare 
various theoretical prescriptions for the first crossing distri- 
bution against our numerical results. Additionally, we wish 
to emphasize the point that the sharp-fc results for the first 
crossing distribution are also insensitive to the non-Gaussian 
process. We therefore show the results for sharp-fc separately 
in this subsection, with those for Gaussian filtered walks in 
the next subsection. 

Figure [2] shows the results for — 0.3, with the red 
(upper) and blue (lower) histograms showing the first cross- 
ing distributions of the constant barrier S = 8 C , for walks 
with positive and negative ei, respectively. The solid his- 
tograms are for model A and the dotted for model B. The 
histograms clearly demonstrate that, within the numerical 
accuracy we reach, the sharp-fc first crossing distribution is 
insensitive to the details of the non-Gaussian process. 

The value |ei| = 0.3 corresponds to a rather large non- 
Gaussianity, since e.g. in the local model one has ei ~ 
3 x 10~ 4 /nl, so that in our case |/nl| ~ 1000. While 
we made this choice to visually enhance the effect of non- 
Gaussianity in the first crossing distribution, it also raises 
an interesting question regarding comparison with theory. 
As discussed elsewhere (see e.g. D'Amico et al. 2011a and 
PGH11), most theoretical analyses rely on a perturbative 
expansion in parameter combinations such as t\v, and the 
large value |ei| = 0.3 means that these analyses break down 
around v < 3. The model presented in PGH11, on the other 
hand, was specifically constructed to remain well-defined in 



Figure 3. Same as Figure[2] but for |ei| = 0.1, with histograms 
shown only for model A. For this case, we also show the theoretical 
predictions of D'Amico et al. (2011a) as the dashed curves. 

the tail v 3> 1 for ei > 0, and is therefore a natural choice 
for comparison. 

The solid curves show the theory predictions from 
PGH11 (i.e., 2 times the expression in Eq. (0) for posi- 
tive (red, upper) and negative (blue, lower) ei , respectively. 
We see that this prescription performs remarkably well when 
ei > 0, while it does not do well for ei < 0. The latter is not 
surprising, given that the theoretical derivation of PGH11 
strictly works only for e\ > 0, with the ei < prescrip- 
tion being ad hoc. The agreement for ei > 0, on the other 
hand, highlights once more the insensitivity of the first cross- 
ing distribution to the details of the non-Gaussianity. This 
is because the calculation of PGH11 assumed the relations 
£n = (ei) n , whereas in our models A and B one can checlQ 
e.g. that £2 ei- 

The dotted curves in Figure [2] show the prediction of 
MVJ00 (i.e., 2 times Eq. ©). Given the large value of |ei|, 
a priori one does not expect this prescription to be a good 
description of the walks, for either sign of t\. This is quite 
clearly the case for ei > 0, with the dotted curve overesti- 
mating the numerical answer by a factor of 4 or more. Re- 
markably, though, the MVJ00 prescription works quite well 

4 Model A has a formal pathology, in that £2 logarithmically di- 
verges; in terms of a cutoff at some s = So, the relation at leading 
order in e\ is £2 = e?(26/18+ln(S/So)/6). Our numerical solution 
has a natural cutoff So = As, although in principle one can regu- 
late the divergence more carefully. For our walks, the logarithmic 
term at most contributes at the same order as the constant 26/18. 
We chose not to be more elaborate with the regularization, since 
model B shows results identical to model A, without having any 
such pathology; in model B, at leading order we find e2 = (4/3)e?. 
Both models have 62 ^ ef, and demonstrate the robustness of the 
first crossing distribution against changes in the exact form of 
this relation. 
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Figure 4. First crossing distribution for Gaussian filtered walks having non-Gaussian initial conditions specified by model A (histograms) 
and with Pg at", with n = —1.2 (left panel) and n = — 2 (right panel). The red (upper) and blue (lower) histograms show results for 
e\ = 0.3 and ei = —0.3, respectively. For comparison, the black (middle) histograms show the first crossing distribution for corresponding 
walks with Gaussian initial conditions. The smooth curves are the same in each panel. The black solid (middle) curve is the Gaussian 
prediction Eq. The other solid curves are the non-Gaussian prediction Eq. 10, with the red (upper) and blue (lower) corresponding 
to e± = 0.3 and e± = —0.3, respectively. The dotted curves similarly show the prediction Eq. 0. See text for a discussion. 



for ei < 0, although it slightly overestimates the numerical 
answer. 

Figure [3] has the same format as Figure [21 and shows 
the sharp-A; results for a smaller value [ei | = 0.1. This time, 
in addition to PGH11 and MVJ00, we also show the cal- 
culation of D'Amico et al. (2011a) (dashed lines). We see 
that all three prescriptions perform reasonably well now, al- 
though once again the MVJ00 prescription seems to work 
better than the others at negative ei. It is not clear to us 
whether this is indicating something deeper, or is simply a 
coincidence. 



3.4 Non-Gaussian walks: Gaussian filtered case 

We now turn to the more interesting case of filtered walks. 
Figure [4] has a similar format as the previous figures, and 
shows our results for Gaussian filtered model A walks, for 
the first crossing of the constant barrier S — 5 C . The red 
(upper) and blue (lower) histograms are for ei = ±0.3 as 
before, while the black (middle) histogram shows the an- 
swer for Gaussian initial conditions, with the corresponding 
analytical prediction from Eq. (3J. The two panels are for 
n — —1.2 and —2, and clearly the corresponding histograms 
do not show any significant difference from each other, high- 
lighting the point we made earlier regarding the irrelevance 
of the form of the power spectrum for the completely corre- 
lated tail. We expect the same to remain true for ACDM as 
well. 

The theory curves in each panel are the same; the solid 
curves show Eq. © (red, upper for ei = 0.3 and blue, 
lower for ei = —0.3), while the dotted curves similarly show 



Eq. ([8]). Notice that these curves describe the histograms in 
each panel as well (or poorly) as the corresponding curves 
in Figure [2] described the sharp-fc non-Gaussian walks. This 
is a clear demonstration that the factor 1/2 argument works 
extremely well in the tail of the distribution, independently 
of choice of power spectrum. 

The histograms in Figure [5] show the results for Gaus- 
sian filtered model B walks with n = —1.2. For the theory 
curves, we only display Eq. J7J for ei > (red, solid) and 
Eq. (SJ for e\ < (blue, dotted). Additionally, the dashed 
lines show the prediction in Eq. (|13|) of single-point statis- 
tics for model B. The dashed lines compare extremely well 
with the histograms, and confirm our argument that single- 
point statistics are sufficient to describe the first crossing 
distribution of filtered walks. The good agreement between 
the dashed lines and the other theory curves indicates that 
the factor 1/2 argument is also working well, and is more- 
over independent of the type of non-Gaussianity. To make 
this point explicitly, Figure [6] shows the ratio of the (nu- 
merical) first crossing distributions for the Gaussian filtered 
and sharp-fc walks. We see that the factor 1/2 argument 
works well within our numerical accuracy. There is an order 
~ 10% discrepancy at larger s for the case of n — —2, model 
A (middle panel), indicating the limit until which the walks 
can be assumed to be completely correlated. 



4 DISCUSSION 

An essential ingredient in the application of the excursion set 
approach to study structure formation from non-Gaussian 
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Figure 5. Same as Figure[4] but for non-Gaussian walks defined 
by model B, with n = —1.2. For the theory curves, we only display 
the ei > prediction for Eq. J7J (solid red, upper) and the t\ < 
prediction for Eq. JS} (dotted blue, lower). Additionally, the 
dashed curves show the prediction of single-point statistics for 
model B in Eq. (1 1 3 I t for positive (red, upper) and negative (blue, 
lower) ex, respectively, with |ei| = 0.3. See text for a discussion. 



initial conditions, is the ability to calculate the first cross- 
ing distribution /ng(s) of a chosen barrier, for the random 
walk performed by the non-Gaussian initial density field 8(s) 
smoothed on a scale corresponding to variance s. We pointed 
out that, while there exist several prescriptions in the litera- 
ture to calculate /ng(s) theoretically, as yet there has been 
no systematic test of these prescriptions using numerically 
generated first crossing distributions. We have performed 
such tests in this work, and in the process have gained some 
insight into the structure of /ng(s). We summarize our re- 
sults below: 

• By simulating two models of non-Gaussianity (NG) 
(Eqs. 1111 and I12[) , we argued that for the case of sharp-fc 
filtering, /ng(s) is insensitive to the exact details of the 
non-Gaussian process. We also showed that for small values 
of reduced skewncss eij, three different theoretical prescrip- 
tions (PGH11, MVJ00 and D'Amico et al. 2011a) work rea- 
sonably well (Figure [3}, while for larger values of |ex|, the 
PGH11 prescription works well for e± > while the MVJ00 
prescription does better for e\ < (Figure [5}. 

• For the more interesting case of filtered walks, we ar- 
gued that the small s tail of the first crossing distribution 
is essentially determined by the single-point statistics of the 
smoothed non-Gaussian density field. For one of our non- 
Gaussian models (model B) which admits a closed form an- 
alytic expression for the single-point PDF, we also explicitly 
verified this numerically (Figure [5}. 

• We further argued, and demonstrated numerically (Fig- 
ures U and [SJ) , that the first crossing distribution for filtered 
non-Gaussian walks in the limit of small s, is simply a factor 
1/2 times the corresponding distribution for sharp-fc walks. 




Figure 6. Behaviour of the ratio /filtered //sharp-fc of Gaussian 
filtered and sharp-fc first crossing distributions, relative to the 
predicted value of 1/2. The three panels show various combina- 
tions of power spectrum Ps(k) oc k n and non-Gaussian process 
as discussed in the text. The black crosses, red squares and blue 
triangles show the results (with Poisson errors) for ei = 0, 0.3 
and —0.3, respectively. 



(For the numerical results, we used Gaussian filtering, but 
the results are not expected to change for TopHat filtering.) 
Since walks with Gaussian initial conditions show the same 
effect, this implies that the ratio /NG//gauss is independent 
of choice of filter, providing a theoretical motivation for the 
practice of prescribing such a ratio when dealing with non- 
Gaussian mass functions. 

• We demonstrated that the previous result is remarkably 
robust against changing not only the power spectrum (which 
might have been expected based on results for Gaussian ini- 
tial conditions), but also the details of the non-Gaussian 
process. While this probably means that halo abundances 
on their own will never be very sensitive to differences in 
types of NG, it does also reassure us that theoretical cal- 
culations can safely make assumptions regarding the non- 
Gaussian process (e.g., the PGH11 assumption of e_, = (ei) J , 
or simply the single-point statistics of model B in Eq. (|13[l ). 
without significantly affecting the final answer. 

These considerations imply that studying primordial NG 
in the framework of the excursion set approach is actu- 
ally much simpler than sharp-fc analyses seem to suggest: 
at small s, the first crossing distributions of filtered walks 
(which are expected to be more realistic than sharp-fc walks) 
are determined entirely by the single-scale statistics of the 
initial density field, with no complications due to the cloud- 
in-cloud problem. To the extent that first crossing distri- 
butions are sufficient to study halo abundances, we would 
argue that, for large masses, it is both simpler and more ac- 
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curate to study primordial NG using a Press-Schechter-like 
approach rather than the full machinery of excursion sets. 

However, as mentioned earlier, equating a first crossing 
distribution to a halo mass function is only an ansatz. All 
our results are subject to the same caveats as in the case of 
Gaussian initial conditions, which were discussed by PLS11. 
We will therefore not repeat them here, except to note that 
the issue of stochasticity in the barrier height as a function 
of scale (or simply a stochasticity in the constant value of 
So) will affect the functional form of the ratio /NG//ga USS - 
E.g., as PGH11 noted, simply prescribing a ratio with some 
chosen value of 5 C is not sufficient, since the behaviour of the 
mass function in the tail will then depend on the choice of 
Gaussian mass function. In particular, not all combinations 
of Gaussian mass function and S c values in the ratio will 
lead to well-behaved non-Gaussian mass functions. Address- 
ing this issue, however, will require a better understanding 
of the jump between first crossing distributions and mass 
functions. 

One issue which we have not addressed in this work is the 
question of halo bias in the presence of non-Gaussianities. 
Paranjape & Sheth (2011) recently discussed halo bias in 
the excursion set framework with correlated steps, for Gaus- 
sian initial conditions. Their results built on a simple scal- 
ing ansatz proposed by PLS11 for writing the conditional 
first crossing distribution, which was inspired by bivariate 
Gaussian statistics and was shown to work rather well when 
compared with Monte Carlo simulations. It will be interest- 
ing to see whether a similarly simple ansatz can be devel- 
oped in the non-Gaussian case, especially since, as we have 
seen, the effect of the filter at least is exactly the same here 
as in the Gaussian case. Previous studies of halo bias with 
non-Gaussianities suggest that, in this case, the choice of 
non-Gaussian process should be far more relevant than it 
was for the first crossing distribution. We leave these issues 
to future work. 

Another discussion we did not enter was the problem of 
voids in the excursion set approach (Sheth & van de Wey- 
gaert 2004; Kamionkowski, Verde & Jimenez 2008; Lam, 
Sheth & Desjacques 2009; D'Amico et al. 2011b). For the 
case of Gaussian initial conditions, Paranjape, Lam & Sheth 
(2011b) recently argued for a formulation of the problem in 
terms of two barriers, one of which is constant and negative 
and the other linearly decreasing from positive to negative 
values. The first crossing problem requires counting those 
walks which first cross the constant barrier, without hav- 
ing crossed the linear barrier at smaller s. In this case, the 
difference between sharp-fc and filtered walks is dramatic; 
whereas the sharp-fc case is extremely sensitive to the exact 
details of the falling barrier, the filtered case is very robust, 
and in fact, the required first crossing distribution is very 
well approximated by the completely correlated answer for 
the single constant barrier, Eq. ©. This latter result in 
particular is encouraging for the case of non-Gaussian ini- 
tial conditions; we expect that the void first crossing distri- 
bution in the presence of non-Gaussianities should also be 
well described by the single point statistics of the smoothed 
non-Gaussian density field using a single (negative) constant 
barrier. However, the sensitivity of the sharp-A: result to the 
details of the problem implies that the ratio prescription will 
probably not work well in this case. 

As a final remark, we recall that PLS11 showed, for Gaus- 



sian initial conditions, that the first crossing distribution of 
a constant barrier by completely correlated walks can be 
trivially extended to the case when the barrier is a mono- 
tonic function B(s). Almost all that was needed there, was 
a replacement S c — > B(s) in the survival probability. In the 
non-Gaussian case, since we are only interested in the small 
s, completely correlated tail of the first crossing distribu- 
tion, and since the arguments of PLS11 did not rely on the 
initial conditions being Gaussian, it should be a straightfor- 
ward matter to extend the non-Gaussian constant-barrier 
results to the case of moving barriers. In particular, this 
seems to be a promising way of approaching the problem of 
re-ionization in the presence of non-Gaussianity, following, 
e.g., Furlanetto, Zaldarriaga & Hernquist (2004). 
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APPENDIX A: THE UNEQUAL SCALE 
BEHAVIOR OF THE THREE-POINT 
FUNCTION 

In this Appendix we investigate the behavior of the three- 
point function {51S2S3) where 8j = 5(sj) = 5(s(Rj)), in 
the regime where Ri 3> R2 2> Rs, i.e. when S(Ri) <C 
S(R2) <C S(Ri). We begin by discussing models inspired 
from inflationary physics, in particular the local model, and 
then compare with the simpler models we have considered 
in this work. 

Al The local model 

The three-point function of S(R) can always be written in 
terms of the bispectrum B(ki, fca, ks) of the non-Gaussian 
curvature pertubation £ generated during inflation, as 

<<5i<5 2 <5 3 ) oc J dkiklT{ki)W(kiRi) J Ak2klT{k 2 )W(k 2 R2) 

x f 1+2 dyy 3 T(y)W(yR 3 )B(k 1 ,k2,y), (Al) 
J fc 2 — fei 

where T(k) is the matter transfer function. For a CDM cos- 
mology, T(k) decays like T(k) ~ (l/fc) 2 logfc for k > fc cq , 
where fc oq is the scale corresponding to matter-radiation 
equality (Bardeen et al. 1986). 

With the assumed hierarchy of scales, this expression 
for {616283 } is never very sensitive to -R3. Indeed, because 
of the filters W(kiRi) and W^fe-fi^), the integrand is sup- 
pressed much before k\ + fc 2 becomes of order I/R3, and 
one can thus assume y <C I/R3, for which W{yR3) ~ 1. 
Moreover, since I/R2 3> 1/Ri, the fc 2 integral has a much 
larger range than the ki one. For k 2 2> ki the width of 
the y integral becomes small, and its value can be approxi- 
mated by the value of the integrand at fc 2 times the width, 
i.e. 2ftife|T(fe 2 )B(fei, k 2 , k2)W(k 2 R3)- 

For non-Gaussianity of the local type, the bispectrum 
gives B(ki,ka,fa) oc 2P c (fci)P c (fc 2 ) + P 2 (fc 2 ), where P ( (k) is 
the primordial power spectrum of curvature perturbations. 
In the fc 2 ^> ki limit, the integrand above thus reduces to 
the sum of the two factorized terms 

4 [ktP(ki)T(k!)W(kiRi)] 
x [k^P(k2)T 2 (k2)W(k2R2)W{k 2 R3)] , (A2) 

and 

2 [kfT{ki)W{kiRi)] 

x [ktP 2 {k2)T 2 {k 2 )W(k2R2)W{k2R3)\ ■ (A3) 

For large fe 2 (and assuming a scale invariant power spec- 
trum), the first term scales like (1/fe) log 2 (fc 2 ), which would 
diverge if it were not for the filter. The integral over fc 2 is 
thus dominated by the value around I/R2, and is largely in- 
sensitive to the exact behavior of the integrand for fc 2 < k%. 
The second term instead scales like (l/fcf) l°g 2 (fe), which is 
convergent and subleading. The integral in (|A1|) is therefore 



well approximated by the factorized integral of (|A2|) over k\ 
and fc 2 , the second of which is simply the two-point function 

{8283) oc j Ak 2 ktP(k2)T 2 (k2)W{k2R2)W{k2R3) ■ (A4) 

Also, the integral over ki is dominated by values around 
1/Ri, where its integrand (at scales smaller than the equal- 
ity scale) behaves like (1/fei) log(fci). Confronting it with the 
variance S(Ri), which in the same regime scales like the in- 
tegral of (1/fci) log 2 (fci), we conclude that the scaling of the 
three-point function is 

( 8(S 1 )6(S2)S(S3) ) ~ VST < 5{S 2 )8(S 3 ) ) , (A5) 
which reduces to VS1S2 with a sharp-fc filter. 

A2 Models A and B 

Turning now to the non-Gaussian models used in this paper, 
the three-point function of model A defined in Eq. © reads 
at linear order in ft 

/ r r c \ g(gl) [ S \, (S2S 1 )(S3S(S)) 

{S1S2S3 } = \ / ds- ' + 5 perms.; (A6) 

3 Jo 7(«) 

from this expression, setting Si = S2 = S3 — S one gets the 
skewness 

( 5 3 (S) ) = 2/3(S)S / 5 ds <^M( g )) (A7 ) 
Jo 7(«) 

and setting ( 5 :i ) — eiS 3 ' 2 , where ei is the reduced third mo- 
ment of the model we want to reproduce, one immediately 
recovers Eq. (|10|) . 

The equal scale behavior of the three-point function can 
be matched exactly by fixing ft, for an arbitrary function 
7(5). We now want to use the residual freedom to reproduce 
(at least approximately) the unequal scale behavior. In the 
sharp-A: case one has (<5j<5j ) = mm(Si,Sj) and (5i5(s)) = 
•&(Si — s); setting y(s) = 4-^/s simply yields 

(W 3 ) = ^5 1 3/2 + ^±^(S 1 ^ + 5 2 ^) (A8) 

with Pi — ei(Si). We notice that, as long as ei is nearly 
constant, this result is nearly independent of S3. Secondly, 
in the 5*i <C S2 limit the dominant term is the one with 
and the three-point function becomes 

<<5i«5 2 «5 3 > =i jVS\S 2 , (A9) 

which matches the approximate scaling of the local model. 

For a generic filter the computations are slightly less 
immediate, but one can check that the very same argu- 
ments that lead to the previous formula go through, with 
S2 replaced by (<5 2 <53), just like what happens for the lo- 
cal model. This matching could in principle be made more 
accurate by choosing 7 so that the exact scaling with Si of 
the ki integral in Eq. (|A2|I is reproduced. However, since our 
simulations show that the fine details of the unequal scale 
behavior of ( 5i<5 2 <53 ) are not very important, for the level of 
accuracy requested in this paper we are satisfied with the 
simple expression for y(s). 

Finally, the three-point function of the non-Gaussian 
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model B defined in Eq. (|12[) . for a generic filter, reads 
<<5iW = ~(SiS 2 S 3 ) 1/2 



/3iri2ri3 + /32r2ir 2 3 + /?3»'31»'32 



+ g(/3l/32^3)'"12r23T'31 



(A10) 



{5 i S j )/i/S~Sj, and ft = /3(s 8 ). For 



where we defined r; 
the sharp-fc filter, at leading order in /3, this becomes 

2 



,9i5' 1 3/2 + /?25'iVS2"+/33 5l ^ ; 



(All) 



Again, setting Si = S2 = S3 = S and matching against 
(5 3 ) = ei5' 3/2 yields = ei/2. In the limit Si < S 2 < 

S3 the dominant term is the second one: we have therefore 



' Si 82S3 ) 



£l 



S1VS2 



(A12) 



which does not match the scaling of the local model in the 
hierarchical regime. Similar results follow for generic filters 
by using Eq. (|A1D|I . 



APPENDIX B: PDF FOR MODEL B 

The non-Gaussian PDF for v = 8/VS in Model B can be 
computed from the Gaussian PDF for Sg as 

f e~ x2/2 ( x 2 -l\ 
p(u) = dx — y^^Su \v — x — e — - — , (Bl) 



2tt 

where e = ei. One now applies the relation 8o[f(x)] — 
J^. So{x — Xi) /\f' (xi)\ (the XiS being the roots of the equa- 
tion fix) = 0) for the change of variables in the Dirac delta 
function. Here f(x) = is just a quadratic equation with 
the two roots 



3 

x± = - 
e 



(B2) 



while its Jacobian gives 



the PDF therefore reads 



-/2 + e -*±/2 



P(v) = — . (B4) 

^F x /l + (2e/3)(i/ + e/6) 

For small e one has x+ ~ z/ and a;_ ~ — 6/e; therefore in this 

_ 2 /9 — 2 /9 

limit e -' is exponentially suppressed, and from e +' 
one correctly recovers the Gaussian result. For generic val- 
ues of e, plugging in the above expression the roots x+ and 
x- from Eq. (|B2[) gives the result quoted in Eq. (Q3J). The 
first crossing distribution for walks with completely corre- 
lated steps follows from constructing the survival probabil- 
ity j Sa l^ &vp(v) and then differentiating with respect to s, 
which leads to the expression in Eq. (|13[) . 



